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ABSTRACT 

We present a detailed single pulse study of PSR B1055-52 based on observations at 
the Parkes radio telescope. The radio emission is found to have a complex modulation 
dominated by a periodicity of ^ 20 times its rotational period P (0.197s), whose phase 
and strength depends on pulse longitude. This periodicity exhibits a phase-locked de- 
lay of about 2.5P between the main pulse (MP) and interpulse (IP), presumed to be 
the opposite poles of the pulsar. This delay corresponds to a light travel distance of 
many times the light cylinder radius. More complex modulations are found within the 
MP on timescales down to about 9P, and both these and the principal modulation vary 
strongly across the (at least) 7 components which the MP and IP exhibit. The nature 
of the single pulse emission, which ranges from smooth and longitudinally extended 
to 'spiky', is also component-dependent. Despite these disparities, the total pulse in- 
tensity distributions at the MP and IP are virtually identical in shape, suggesting a 
common emission mechanism. In an attempt to account for the complex modulations 
we examine a number of physical models, both intrinsic (which presuppose the pulsar 
to be an isolated neutron star) and extrinsic (appealing to the presence of circumstel- 
lar material to modulate the emission). Significant objections can be made to each 
model, so this pulsar's behaviour patterns remain a crucial challenge to theorists. 
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1 INTRODUCTION 

PSR B1055-52 is a middle-aged pulsar (characteristic age 
about 0.5 Myr) known to emit pulsed emission from the 
radio to the 7-ray band. At radio frequencies it exhibits a 
complex integrated profile with separate main pulse (MP) 
and interpulse (IP) contributions, separated by roughly 



180° in pulse long 


itude and arj 


^ued to come from separate 
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: WanE et al.l I2OO6I"). In an earlier paoer 



pattern and polarization of the integrated profile were ex- 
ploited to develop a self-consistent geometric model within 
which the MP and IP emissions come from regions as- 
sociated with the opposing poles of the pulsar's magnetic 
axis, which is inclined to the rotation axis by about 75°. 
PSR B105 5-52 has rec ently been confirmed as a pulsed 7- 
ray source (IBecker fc Tr iimper 1997; Thompson et al., 1999i : 



lAbdo et al.fl2010l ). The pulsed high energy components are 



then seen as emanating from higher altitudes in an outer 
gap. 
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An unexpected conclusion of Paper I was that, while the 
IP emission comes from open field lines and its profile spans 
the appropriate range of pulse longitudes, our oblique sight- 
line towards the pulsar requires that the more complex MP 
emission cannot originate above a conventional polar cap. 
In fact, assuming the linear polari zation can be interpreted 
with the rotating vector model (|Radhakrishnan "fc" Cookd 
Il969h and a strict dipolar field throughout the magneto- 
sphere, the MP appears to be generated on field lines nor- 
mally associated with the pulsar's closed corotating zone, 
suggesting that the true corotating zone is considerably 
smaller. Furthermore, the deduced location of the MP is 
asymmetrically placed wi thin this zone. T he recent im- 
proved 7-ray pulse profile (|Abdo et al.ll2010l ) spans the MP 
in its trailing half and an outer gap model is shown by these 
authors to be consistent with the geometry deduced from 
Paper I. 

Our purpose here is to link our single-pulse radio data 
of PSR B1055-52 with the two-pole viewing geometry elu- 
cidated in Paper I. A pioneering study of this pulsa r (at 645 
MHz) was first published 20 years ago l|Biggs|[l990l ). identi- 
fying the principal pulse profile components and attempt- 
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Figure 1. The pulse profile of PSR B1055-52 at a wavelength of 20 cm showing total intensity (solid line), linear polarization (dashed 
line) and circular polarization (dotted line) reproduced from Paper I. The component numbering convention used in this paper is 
indicated. 



ing to resolve the basic pulse modulations. Evidence was 
found for periodic modulations, delayed correlation between 
the MP and the IP emission and for long-term changes in 
the profile at his observing frequency. Our intention here 
is to exploit the improved signal-to-noise ratio of our data 
to clarify and refine these results, albeit at a significantly 
higher observing frequency. We will demonstrate that each 
component has its own characteristic single-pulse behaviour, 
while exhibiting a subtle level of MP-IP interaction. As far 
as possible, we will attempt to indicate how these results, 
combined with the geometric constraints derived in Paper I, 
give insight into the physics of the pulsar system. 

For the reader's reference we reproduce the image of 
the integrated profile at 20 cm wavelength presented in Pa- 
per I (Fig. [U where we were able to identify an additional 
component in the MP leading to a modified version of the 
component nomenclature introduced by Biggs. 

Some years ago an exercise similar to the present 
investigation was unde rtaken with the pulsar B1702-19 
l|Wehevrede et al.ll200'iD . which, like PSR B 1055-52, has rel- 
atively strong IP emission. In that pulsar the MP and IP 
were much narrower than those of PSR 1055-52 and the 
MP-IP separation was much closer to 180° . This was shown 
to result from the rotation of a highly inclined dipole field. 
Single pulse data modulations of ~ IIP were found in both 
the MP and IP emission which, crucially, exhibited an inter- 
pole time delay of about half a rotation period. Communica- 
tion between the polar regions is a fundamental clue to the 
underlying mechanism whereby the the pulsar emits. One 
may ask whether the communication is carried through the 
magnetosphere, the body of the neutron star or even via 
some external material. Our purpose here is to investigate 
whether PSR B 1055-52 can help to answer some of these 
questions. 

In the next section we present the record of our obser- 
vations, followed by (in Section [3]) a detailed account of the 
components' individual modulation. In Section|4]we discuss 



the issues arising from our results and consider alternative 
models. Finally in Section [S] we summarise our conclusions. 



2 OBSERVATIONS 

The observations of PSR B 1055-52 presented in this pa- 
per were all recorded at the Parkes radio telescope in Aus- 
tralia using different backends and receiver systems. The 
main conclusions of this paper are based on a two hour long 
observation recorded on 2007 April 23. All measurements 
presented in this paper are based on this observation, unless 
stated otherwise. The observation made use of the H-OH 
receiver centred at a frequency of 1390 MHz with a band- 
width of 256 MHz. The signal was processed on-line using 
the analogue filterbank, which combined the two polariza- 
tion channels of the linear feeds of the receiver to give the 
total intensity (Stokes I) sampled once every 50/xs for 512 
separate 0.5 MHz wide frequency channels. The data were 
recorded on tape for further off-line processing. 

A second observation presented in this paper is a 1.9 
minute long observation which was recorded on 2008 July 25 
during observations for the Fermi timing program at Parkes 
(Weltevrede et al. 2010). During this observation the pulsar 
was brighter than usual, most likely because of interstellar 
scintillation. The data were recorded in a very similar way 
compared to the long observation, except that the data were 
sampled at SO/is and that the centre beam of the multibeam 
receiver was used at a centre frequency of 1369 MHz and a 
bandwidth of 256 MHz. 

The third observation we made use of is an hour long 
observation recorded on 2008 October 24. As for the short 
observation the receiver used was the centre beam of the 
multibeam receiver at a centre frequency of 1369 MHz and 
a bandwidth of 400 MHz. The dat a were recorded using the 
BPSR backend (|Keith et al.ll2010l ) with a sampling time of 
64/is and a spectral resolution of 400 channels. The observa- 
tion was shortly interrupted to allow the feeds of the receiver 
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to be rotated by 90° . The two halves were appended off-hne 
to make a single one hour long observation. Unfortunately 
the pulsar was relatively weak during this observation. 

The first step in the off-line data reduction involved the 
masking of frequency channels affected by radio frequency 
interference (RFI). The remaining frequency channels were 
added together accounting for the dispersion delay caused 
by the interstellar medium using the SIGPROC software 
packag43- Using the known pulse period the resulting dedis- 
persed time series were folded following the procedure out- 
lined in lWeltevrede et al.l (|2006t ). The data are now arranged 
in a two dimensional array (pulse-stack) where the rota- 
tional phase (pulse longitude) of the star is one axis and the 
pulse number the other. 

Although narrow band RFI was removed from the data, 
there was still some impulsive broad-band RFI affecting the 
2 hour observation. The data were therefore inspected by 
eye and 2% of the pulses in the pulse-stack were flagged and 
ignored during further analysis. 

The next step in the off-line data reduction attempts to 
correct what we will refer to as "baseline artifacts". These 
artifacts arise when relatively strong signals are recorded us- 
ing the analogue filterbank which records the data with 1-bit 
sampling. The result is that the off-pulse intensity (baseline) 
is not constant throughout the observation, but is decreased 
both before and after a strong signal is recordecQ The base- 
line artifact is not severe in the two hour long observation, 
but it is much stronger in the 1.9 minute observation dur- 
ing which the pulsar appeared significantly brighter and for 
which the sampling rate was slightly lower. 

As a result of the baseline artifact adjacent pulse longi- 
tude bins are not entirely independent of each other, which is 
undesirable. The baseline artifact can be partially removec|f| 
by convolving the recorded (dedispersed) signal with two ex- 
ponential functions, one removing the decrease in the noise 
level after a strong signal and one removing the noise level 
before a strong signal. After this procedure the off-pulse 
baseline no longer showed the characteristic dips directly 
before and after the pulse. Finally a constant baseline was 
subtracted from each pulse to make the average intensity in 
the off-pulse region zero. 

Although we do not expect the adjacent pulse longitude 
bins to be completely independent of each other after the 
correction, at least it gives us a tool to test if the baseline 
artifacts can be expected to influence our results. We find 
no appreciable differences in any of the results presented 
in this paper when using the corrected or uncorrected data. 
This suggests that the correlations that will be discussed are 
intrinsic to the pulsar rather than instrumental. In addition 
the BPSR data, which did not suffer from baseline artifacts, 
are fully consistent with the results presented in this paper. 



^ http:/ /sigproc. sourceforge.net/ 

^ The reason why the noise level before a strong signal in a 
frequency channel can be affected is the combination of the earlier 
arrival of the strong signal at frequency channels corresponding 
to higher sky frequencies in combination with frequency channels 
sharing filters. 

^ It is impossible to remove the baseline artifact perfectly, be- 
cause of its complex time-frequency-intensity behaviour. 
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Figure 2. Left: A pulse-stack of 581 successive pulses of the MP 
observed when the pulsar was exceptionally bright. Right: The 
pulse-stack of the same pulses as in the left panel, but now centred 
at the IP. The ~ 20P modulation pattern is visible in the central 
components of the MP and across the IP, as is the delay of ^ 2.5P 
by which the pattern observed in the MP lags the IP. 



3 SINGLE PULSE MODULATION 
3.1 Pulse-stack 

Before embarking on a detailed analysis of the statistical 
properties of the emission of PSR B1055-52, it is good to 
first get a basic understanding of its behaviour by inspect- 
ing some of its individual pulses by eye. The pulse-stacks of 
the MP and IP of PSR B1055-52 are plotted in Fig.[lwith 
the pulse longitude scale aligned with that of Fig. [T] by plac- 
ing the profile peak at pulse longitude 90°. As one can see, 
the pulse intensities are highly variable from pulse to pulse. 
The most obvious modulation pattern is the roughly peri- 
odic brightening which happens on a timescale of ~ 20P. 
This modulation is especially visible in the IP and the third 
component of the MP (at ~ 100° pulse longitude). 

Periodicities in pulse i ntensity are a common phe- 
nomenon in radio pulsars (e.g. lWeltevrede et al. 2006, 200"^) , 
often associated with drifting subpulses (jPrake fc CraftI 
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Il968l ). although in the case of PSR B1055-52 the periodicity 
arises from intensity modulations rather than phase modu- 
lations. We will show in Sect. 13.71 that there exists a stable 
phase delay between the modulation patterns observed in 
the MP and IP, a feature which in fact is detectable by eye 
in Fig. 

The various components of the pulse profile have com- 
pletely difi'erent characteristics in their single pulse be- 
haviour. It is particularly striking that the pulses of the 
MP have a very difi'erent appearance to those of the IP. 
The pulses of the MP show much sharper structures and 
can have widths much smaller than the profile component 
widths. In contrast the individual pulses of the IP are much 
more extended or smeared out in pulse longitude. The pulses 
of the MP, especially those of component II (pulse longitude 
~ 90°), can be characterised as very spiky, i.e. the pulse- 
stack at the location of component II is dotted with nar- 
row spikes of emission with high peak fiuxes thereby resem- 
bling the emission obser ved, for instance, in PSR B0656-f 14 
l| Weltevrede et al]|2006l ). Furthermore, it appears that not 
every MP component shares the 20P modulation as clearly 
as some of the others. 

Also within the IP the individual components have dif- 
fering single-pulse characteristics. Although component VI 
(pulse longitude ~ 248°) participates in the modulation cy- 
cle, it never fades as much as the other IP components, 
thereby forming a vertical band in the pulse-stack. In addi- 
tion, component V (pulse longitude ~ 240°) and VII (pulse 
longitude ~ 260°) while also subject to modulation, do not 
completely disappear. The middle component (VI) appears 
to form some sort of boundary in the sense that sometimes 
the leading side of the IP profile participates much more 
strongly in the modulation cycle compared to the trailing 
side (e.g. around pulse 96 and 245) and sometimes this ap- 
pears to be reversed (e.g. around pulse 85 and 355). Nev- 
ertheless in many cases both components clearly show the 
same cycle. The duty cycle of the ~ 20P modulation pattern 
is relatively small in both the MP and IP with the pattern 
being bright for only a few pulses. 

3.2 Fluctuation spectra 

The time-averaged properties of the periodic subpulse mod- 
ulation can be studied in more detail by calculating fiuc- 
tuation spectra. We will fo llow the methods developed by 
lEdwards fc Stapperj ll2002l) and for more details about the 
analysis we also refer to Weltevrede et all (|2006l ). The first 
type of fluctuation spectrum comput ed is the lon gitude- 
resolved fluctuation spectrum fLRFS: iBackeil [l970l l. This 
fluctuation spectrum is obtained by separating the pulse- 
stack into blocks of 256 pulses long. By calculating Fourier 
transforms along constant pulse longitude columns in the 
pulse-stack one gets information about the periodicities 
present for each pulse longitude bin. By averaging the 
square of the spectra obtained for the consecutive blocks 
of data the LRFS is obtained. These spectra are shown in 
the middle panels of Fig. [3l 

The LRFS confirms and clarifies the features noted by 
eye in the pulse-stack. The LRFS of the IP clearly shows the 
~ 20P periodicity, which corresponds to ~ 0.05 cycles per 
period (cpp). The same periodicity is visible in the LRFS 
of the MP (especially in the trailing half of component II 
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Figure 3. Modulation index and fluctuation spectra of PSR 
B1055-52. The left-hand plots are for the MP and the right-hand 
plots for the IP. Upper panel: The peak of the pulse profile (grey 
curve) is normalised to 5 and the longitude resolved modulation 
index is shown as the points with error bars. Middle panel: The 
LRFS. Lower panel: The 2DFS. The data is re-binned to 1000 
pulse longitude bins. 



at 90° and component III at 100°), although there is much 
more additional modulation power at different frequencies 
(or timescales) compared to the IP. This is consistent with 
the ~ 20P periodicity visible most clearly in the pulse-stack 
of the IP. As can be seen in the side panel of the LRFS of 
the MP, the spectral power appears to be a combination of a 
white-noise like plateau combined with a peak correspond- 
ing to the ~ 20P periodicity and a red-noise like spectrum 
corresponding to timescales larger than ~ 5P (correspond- 
ing to frequencies smaller than ~ 0.02 cpp). The different 
profile components show different mixtures of these spectral 
behaviours. 

The two dimensiona l fiuctuation spectrum (2DFS; 
lEdwards fc Stappers|[20o3 ) is a different type of fiuctuation 
spectrum that can reveal additional information about pe- 
riodic subpulse modulation. The vertical frequency axis is 
very similar to the case of the LRFS, which can be expressed 
as P/P3 , where P3 corresponds to the 20P periodicity visible 
in the pulse-stack. The horizontal axis of the 2DFS denotes a 
frequency as well, which can be expressed as P/ P2 , where P2 
corresponds to the horizontal separation of the subpulses. If 
the power were concentrated at positive values of P2 it would 
indicate that subpulses tend to shift (or drift) to later pulse 
longitudes in successive pulses. 

The IP shows the same P3 = 21 ± 2P periodicity (for 
a discussion about erro r bars on fiuctuation frequencies, see 
I Weltevrede et al.ll2006l ) and its 2DFS is perfectly symmet- 
ric about the vertical axis. This indicates that subpulses in 
successive pulses do not drift on average to later or earlier 
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Figure 4. A section of the S2DFS showing the evolution of P3 by 
sUding a 256-pulse long window over 6255 successive IP pulses, 
thereby generating 6000 spectra. 

pulse longitudes, which is consistent with the fact that the 
intensity bands in the pulse-stack of the IP appears to be 
horizontal. The 2DFS of the MP has a much more complex 
shape. As expected from inspection of the LRFS, the 2DFS 
of the MP also shows power corresponding to a wider range 
of periodicities. There is horizontal structure visible in the 
2DFS with power concentrated on frequencies corresponding 
to the component separations within the MP. The asymme- 
try of the 2DFS about the vertical axis indicates that the 
appearance of subpulses in the different components is cor- 
related, although the correlation is complex in nature. The 
possible pulse longitude drift of subpulses participating in 
the P3 ~ 20P cycle will be explored in more detail in Sec- 
tion ET] 

3.3 Stability of the P3 ^ 20P cycle 

The peaks in the LRFS and 2DFS caused by the P3 ~ 20P 
cycle is broad, hence the periodicity does not appear to be 
particularly stable. Although the first signs are not promis- 
ing, we pursued the possibility that a very precise and stable 
clock is responsible for the observed periodicity, although 
modulated by a noise-like process. The presence of such 
a stable signal would have important consequences, for in- 
stance when the p eriodicity is linked to a possible fallback 
debris disks (jCordes fc Shannon. 2008 ). 

The possible presence of a stable periodicity among 
other time-varying signals can be tested by making use of 
the sliding 2DFS (S2DFS: [Servlak et al.| [2o'o9'). The S2DFS 
is a time-resolved version of the 2DFS and is therefore ideal 
to analyse the time-dependence of the P3 ~ 20P cycle. The 
S2DFS is calculated by computing the 2DFS for blocks of 
256 pulses with the start of the block shifting with one pulse 
at the time. In other words, the first 2DFS is computed from 
pulses 1 — 256, the second for pulses 2 — 257, etc. The two di- 
mensional spectra are collapsed over the P2 axis, giving rise 
to a one dimensional spectrum for each block. By stacking 
the resulting spectra for the successive blocks of data one 
obtains the S2DFS (see Fig. [4|. The spectra obtained for 
these "block numbers" are not independent of each other, 
as the time resolution is determined by the size of the slid- 
ing window of 256 pulses. The S2DFS shows a constantly 
changing frequency of the modulation pattern and no fixed 



periodicity is present among the other time-variable fluc- 
tuations. The second conclusion is that the timescale of the 
time- variability of the modulation pattern is short, i.e. a few 
hundred pulses or less. 

Secondly we tested the hypothesis that the observed 
P3 ~ 20P cycle is a very stable periodicity affected by a 
large phase-noise. In other words, could it be possible that 
the bright intensity bands in the pulse-stack appear often 
a few pulses earlier or later then what would be predicted 
by an otherwise strictly periodic cycle? In some ways this 
is analogous to what is assumed for pulsar timing experi- 
ments, where the pulse arrival times are compared to a tim- 
ing model containing well defined periodicities, giving rise to 
timing residuals with relatively small amp litude. By making 
use of the timing software developed by IWeltevrede et ahl 
we tried to model the measured pulse numbers for 
which we observe a bright intensity band in the IP with a 
timing model being a fixed frequency. This frequency was 
initially set to the frequency of the peak in the LRFS ob- 
tained for the same stretch of data of a few hundred pulses 
long which was used to measure the pulse numbers corre- 
sponding to the appearance of the intensity bands. It turns 
out to be impossible to get phase connection of the measured 
arrival times of the intensity bands, even after changing the 
frequency in the timing model. This confirms the suspicion 
that the periodicity of the modulation pattern is unstable, 
rather being a fixed value modulated by phase noise with an 
amplitude less than 20P. 

3.4 Modulation index 

The upper panels of Fig. [3] show, besides the pulse pro- 
file (grey curve), the longitude resolved modulation index 
of the MP and IP. The modulation index is a normalised 
measure of the variability of the pulse intensities at a given 
pulse longitude and is computed from the LRFS (see e.g. 
lEdwards fc Stapperdliool : IWeltevrede et al.|[200^ '). The er- 
ror bars on the modulation index was derived by boot- 
strapping the data (pulse-stack). For each iteration in the 
bootstrap method an equivalent pulse-stack was generated 
by adding random noise samples with an root-mean-square 
(rms) equal to that measured in the off-pulse region. The 
standard deviation of the distribution of modulation indices 
obtained in this way is set to be the (1-cr) error bar. Sim- 
ulations with artificial pulse-stacks indicate that the error 
bars derived in this way are more reliable than the ana- 
lytically estimated erro r ba rs as computed for instance by 
Edwards_ &_Stappers ( 20031 ) and IWeltevrede et al.l (|2006l '). 

Typically one measures modulation indices of around 
0.5 (at observing fre quencies around 1400 MHz; e.g. 
IWeltevrede et akl I2OO6I) at the centre of components with 
the modulation index fiaring up at the component edges. A 
similar flaring up at the proflle edges is seen for this pulsar. 
The most variable pulse intensities (i.e. highest modulation 
index) are seen between components II and III of the MP. 
The value of the modulation index in this region of the MP 
is as large as that o bserved for PSR B0656+14 at a fre- 
quency of 327 MHz (|Weltevrede et all 2006). However one 
cannot direct compare these values because the high mod- 
ulation index of PSR B1055-52 is observed between profile 
components, while that of PSR B0656-I-14 is observed at 
the profile centre (which is what makes the high modulation 
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Figure 5. The pulse profile variability of 2000P long sub integra- 
tions of the MP (left) and IP (right) is quantified by the longitude 
resolved standard deviation (open circles) and the modulation in- 
dex (points with error bars). The peak of the average pulse profile 
(solid line) is normalised to a value of 0.15. 



index of PSR B0656-fl4 unusual). The modulation index 
derived from the BPSR data is similar, suggesting that the 
high modulation index is not caused by the baseline artifact. 

The modulation index profile of the IP is on average 
lower than that of the MP. This is not surprising given the 
much more spiky appearance of the pulse-stack of the MP. 
The modulation index profile shows a dip at the location of 
component VI, which is caused by the presence of steady 
emission in that component (as is seen as the vertical band 
in the pulse-stack in Fig. [2|. 



3.5 Profile variability 

iBiggd jiggd ) reports that the pulse profile of PSR B1055-52 
is variable at relatively long timescales in the sense that the 
peak amplitudes of the different components fiuctuate when 
integrations of 2000P are taken. In order to quantify this 
effect we averaged pulses in blocks of 2000 pulses. Using the 
methods outlined in the previous subsection we calculated 
the longitude resolved modulation index and standard devi- 
ation (after rebinning the data to 200 pulse longitude bins). 
As can be seen in Fig. [H the highest variability (the open 
circles indicating the standard deviation of the intensities at 
a given pulse longitude) is observed in component II (pulse 
longitude ~ 90°). However, after normalisation with respect 
to the average intensity observed at these longitudes (i.e. the 
modulation index; the points with error bars) one can see 
this variability is not exceptional compared to other com- 
ponents in the MP. The emission of the IP is more stable 
resulting in a low modulation index for the profile variabil- 
ity which is comparable to that observed in component IV 
in the MP (pulse longiture ~ 110°). This is consistent with 
the more spiky nature of the emission observed for the MP, 
hence longer stabilisation ti mescales can be expecte d as is 
observed fo r PSR B065 6-H4 (| Weltevrede et al.ll2006l ). 

Unlike [Siggsl (|l990l ). we find that component II h as the 
larges t amplitude throughout the observation, while iBiggd 
1 1990 ) observed component III to have in general the largest 



amplitud«0. We have found no exception to this rule in 67 
observations of this pulsar betwee n 2007 and 2011 for th e 
Fermi timing program at Parkes (| Weltevrede et aT]|2010l ). 
It must be noted that our ob serving wave length (around 20 
cm) is different from that in lBiggsr(|l99(j ) (around 50 cm), 
but the profile evolution bet ween these frequencies is no t 
very strong (see e.g. Fig. 2 in I Weltevrede &: Wright! [20091 ). 
We attribute these differences in the profile to instrumen- 
tal effects in the observations of Biggs (1990 ), which had a 
much lower time resolut ion. Later Parke s radio observations 
as show n by for instance l Wu et al.l (|l99ll ): lvan Ommen et all 
(199?!); [Thompson et al.l (|l999l ) are consistent with the ob- 
servations presented in this paper. 



3.6 Pulse-energy distribution 

The fiuctuation spectra have shown us in some detail about 
periodicities in the pulse intensity. The modulation index 
tells us, in a very limited way, how variable the pulse in- 
tensities are. More detailed information about the intensity 
variations can be obtained by making a pulse energy distri- 
bution. The pulse energy distribution is calculated by inte- 
grating the intensities as recorded in the pulse-stack over a 
pulse longitude range. This has been done for the full pulse 
longitude range of the MP and IP separately, which gives a 
pulse energy Ei for each stellar rotation i. Given that the 
observations were not fiux calibrated, only the values Ei 
relative to each other are meaningful. We therefore choose 
to normalise the pulse energies Ei by (Emp), the average 
pulse energy measured for the MP. Closer inspection of the 
data reveals that the fiux of PSR B1055-52 is not constant 
throughout the observation because of interstellar scintilla- 
tion. The pulsar is weak at the start and end of the obser- 
vation, but has a nearly constant fiux for 35 minutes. The 
energy distributions shown in Fig. [S] are therefore based on 
these 35 minutes, which much more accurately refiect the in- 
trinsic pulse energy distribution of the source. The off-pulse 
energy distributions (noise distributions) are shown as well 
in order to quantify the uncertainty in the pulse energy mea- 
surements (dashed grey distributions). The off-pulse energy 
distribution is calculated by integrating the off-pulse region 
over a pulse longitude window with equal width compared 

to the on-pulse window. 

Using the methods explained in IWeltevrede et ahl 

(l2006t ) we have fit the shape of the energy distribution of 
the IP by a log-normal distribution convolved with the ob- 
served noise distribution. The parameters of the lognormal 
distribution 



P{E) = 



(E) 



2tt(tE 



exp 



(1) 



we fit for are found to be = —0.13 and a = 0.51. This fit 
matches the observed energy distribution of the IP very well 
(see centre panel of Fig. |6} . 

The fit for the IP energy distribution is reproduced in 
the left panel of Fig. |S] (without any scaling). As one can 
see the shapes of the energy distribution observed for the 



* In the component name convention of iBiggj l IlQQOl) this cor- 
responds to component I and II rather than components II and 
III. 
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Figure 6. The pulse energy distribution of the MP (left panel) and the IP (central panel). The on-pulse energy distribution is plotted 
as the black histogram and the off-pulse distribution as the grey dashed histogram. The pulse energies are normalised to the average MP 
pulse energy. The lognormal fit for the IP (solid grey histogram) is shown in both the left and central panel. The pulse energy distribution 
of component III is plotted in the right panel, where the pulse energies are normalised to the average pulse energy in that component. 



MP and IP (centre panel of Fig. (Sjl are very similar. This 
close resemblance is surprising given the completely diflerent 
characteristics of the pulses as can be seen in the pulse-stack 
(see Fig. [2]) and given the conclusion in Section [3. 2 1 that the 
fluctuation spectra of especially the MP can be described 
by a mixture of different emission types. The resemblance 
of the MP and IP energy distribution is confirmed in other 
(shorter) observations. 

There are small differences in the shape of the en- 
ergy distribution of the MP and IP. The IP shows a tail 
beyond Ei>4 (Emp) which is not present in the MP and 
a slight excess of MP pulses with energies in the range 
2 (-Emp) ^ Ei<4 (Emp) - In addition the energy distribution 
of the MP is skewed to slightly lower energies at the low- 
energy part compared to that of the IP. Closer inspection 
reveals that the energy distributions of most of the individ- 
ual components are very similar to the shape of the overall 
energy distribution of the MP and IP. The exception is com- 
ponent III, which shows a clear break in its pulse energy dis- 
tribution (see right panel of Fig.[6)l and component I, which 
is too weak to accurately determine its pulse energy distri- 
bution. The different shape of the pulse energy distribution 
of component III is therefore responsible for the slight dif- 
ference in the overall energy distribution of the MP and IP. 

Motivated by the similarity between the pulse energy 
distribution of the MP and IP, we investigated if the pulse 
energies of the MP and IP are correlated at timescales 
slightly longer than the m odulation cycle. Followin g a simi- 
lar procedure employed bv lWeltevrede et al.l (|2007l ) for PSR 
B1702-19, we started by averaging blocks of 20 pulses. This 
ensured that the modulation pattern is no longer present in 
the signal, which was confirmed by calculating the LRFS. 
MP and IP energies were calculated for each block by inte- 
grating over the full pulse longitude range over which emis- 
sion was detectable in the pulse profile. Finally, by cross 
correlating the MP and IP sequence of energies it was pos- 
sible to investigate if the power output of the MP and IP 
are related. However, there was very little evidence for cor- 
relation between the MP and IP power output on timescales 
larger than the modulation cycle, apart from long-term cor- 
relations caused by the fact that both the MP and IP are 
affected simultaneously by interstellar scintillation. 



3.7 Phase-locked modulation 

The fact that the MP and IP both show a similar periodic 
modulation pattern is already a remarkable result in itself. 
We will explore this similarity further by determining any 
phase relations of the ~ 20P modulation pattern observed 
in the different profile components. The technique used to 
obtain this phase relation as function of pulse longitude is 
based on fluctuation spectra. Firstly the (complex) LRFS 
was calculated for 128 pulses long blocks of data. Instead of 
calculating the square of the complex numbers in the spec- 
trum, as is done to obtain the LRFS as shown in Fig. [3l the 
phases were calculated. By selecting the spectral bin that 
showed the ~ 20P modulation pattern strongest (which was 
at 0.046875 cpp), the so-called subpulse phase as functio n 
of pulse longitude was obtained (|Edwards fc Stappersl200^ ) . 
The complex numbers obtained for the different blocks of 
data are added coherently using an iterative scheme designed 
to maximize the correlation between the phase proflles of the 
different blocks in the on-pulse region. This is achieved by 
allowing an arbitrary phase offset, which is constant with 
pulse longitude, between the different blocks of data. The 
coherent addition gives rise to the complex subpulse mod- 
ulation envelope, with the complex phase corresponding to 
the subpulse phase. Error bars on the subpulse phase were 
calculated by applying a bootstrap scheme similar to the 
one employed in Section [3.41 for the modulation index. 

The result of the calculation is shown in the lower panels 
of Fig. [7] The actual value of the subpulse phase is not an in- 
trinsic property of the pulsar, but phase differences between 
different pulse longitudes are. We have therefore arbitrarily 
chosen the subpulse phase to be zero at the centre of the IP. 
The upper panels of Fig.[7]show, in addition to the pulse pro- 
file, the amplitude of the subpulse modulation correspond- 
ing to the amplitude of the complex modulation en v elop«13 
(see lEdwards fc Stapper'3l2002l : lEdwards et al.ll2003l ). Note 
that this amplitude depends on the coherence of the modu- 
lation pattern (i.e. the stability of the longitude dependence 
of the subpulse phase as function of pulse longitude). The 

^ The normalisation of the subpulse amplitude profile is such 
that for a pure sinusoidal signal with a frequency equal to the 
centre of a frequency bin in the complex LRFS the subpulse am- 
plitude profile matches the pulse profile with its peak amplitude 
normalised to 1. 
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Figure 7. The upper panels (MP on the left, IP on the right) show the strength of the 20P modulation pattern as function of pulse 
longitude (black curve) and the pulse profile (grey curve). The lower panels show the corresponding subpulse phase as function of pulse 
longitude with an arbitrary phase offset applied to make the subpulse phase zero at the centre of the IP. The solid line indicates the 
slope expected if the pattern is strictly in phase in the pulsar frame. 



results obtained are identical in the first and second half 
of the data and for the data which was not corrected for 
the baseline artifact. In addition the BPSR data gave iden- 
tical results, although the pulsar was much weaker during 
that observation. Hence we are confident that the features 
of Fig. [T] both at the MP and IP, are stable and intrinsic 
properties of the star rather than instrumental effects. 

The subpulse phase relation with pulse longitude reveals 
a high degree of complexity, especially for the MP. Concen- 
trating on the IP first, we see that the subpulse phase is 
roughly independent of pulse longitude. This means that 
that the modulation cycle seen in different pulse longitude 
bins of the IP are in phase, hence the intensity bands seen in 
the pulse-stack are to a good approximation horizontal (at 
least on average). This lack of phase drift makes it unclear 
if the physical mechanism responsible for the modulation 
pattern is the same as is operating in pulsars with drifting 
subpulses. At the outer edges of the IP profile the modula- 
tion pattern appears to be slightly out of phase such that the 
subpulse phase decreases away from the central parts of the 
profile. This indicates that the modulation pattern is slightly 
lagging at the centre of the IP. The central peak of the IP 
at pulse longitude ~ 248° does reveal a slight deviation in 
subpulse phase as well. 

In the MP, at the longitudes where the 20P modula- 
tion is strongest (pulse longitudes 88° — 105°), the subpulse 
phase appears to be roughly in phase as well, although re- 
vealing a weak decline corresponding to a gradual drift from 
trailing to leading components. The subpulse phase offset 
between the IP and the central parts of the MP is about 
40°, which corresponds to a delay of ~ 2.5P such that the 
pattern in the MP is late compared to that of the IP. It must 
be noted that if the modulation patterns in two pulse lon- 
gitude bins are intrinsically in phase (in the pulsar frame), 
they will be observed out of phase because later pulse lon- 



gitudes are observed later in time, allowing the pattern to 
change during the rotation of the star. The slope in subpulse 
phase expected to be observed if the modulation is intrinsi- 
cally in phase is indicated by the almost horizontal baseline 
in Fig. [7] continued through from the MP to the IP. We can 
conclude that the ~ 20P modulation cycle of the IP is in- 
trinsically out of phase with the modulation of the MP at 
any longitude and that the observed phase delay of the mod- 
ulation patterns corresponds to an intrinsic delay of ~ 3P. 
As will be discussed in Section |4j this delay corresponds to 
an extremely long light travel distance which is much larger 
than the light-cylinder radius, thereby challenging existing 
theories for pulsar emission to explain this. 

The phase delay between the MP and the IP can be 
demonstrated in a different way by correlating the measured 
pulse energies for successive pulses in the MP and IP. The 
cross correlation function, using the regions in the MP and 
IP where the modulation is strongest, is shown in Fig. [8] The 
~ 20P modulation cycle is clearly visible in the correlation 
function, showing that the modulation patterns in the MP 
and IP are phase-locked independent of any jitter in the 
periodicity of the pattern. The phase offset in the patterns 
gives rise to an offset of the peak of the correlation function 
from zero lag. The shift can be determined by comparing 
the offset between the cross correlation curve and the time 
reversed correlation function, which is set to 5P in the figure. 
This offset is equal to twice the offset of the symmetric curve 
from zero lag and therefore confirms the lag of ~ 2.5P of the 
pattern observed in the MP compared to that in the IP. 

The behaviour described so far is accompanied by a 
number of strong diversions. The modulation pattern in the 
trailing peak of the MP (at pulse longitude 110°) appears 
to be out of phase with that in the more central regions by 
~ 80° , corresponding to a delay of almost 5P. More subtle is 
the change in subpulse phase in component II at the location 
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Figure 8. The cross correlation between the pulse energies of 
parts of the MP and IP where the ~ 20P modulation pattern 
is strongest. The peak of the correlation function (solid curve), 
which is normalised to 1, peaks at a negative lag. The dashed 
curve is the same function, but time reversed with a shift of 5P. 
This formally confirms the phase-locked ~ 2.5P lag by which the 
MP follows the IP independent of any jitter in the 20P modula- 
tion. 



of the strongest presence of the ~ 20P modulation (pulse 
longitude ~ 90°). Although the subpulse phase is approx- 
imately longitude stationary in that pulse longitude range, 
it decreases towards pulse longitude 90° , indicating that the 
pattern is observed early at the pulse longitudes where the 
modulation is strongest. Close examination of the MP pro- 
file in Fig. [7] reveals that component II is in fact composed 
of two components which overlap. The trailing half of com- 
ponent II shows the 20P modulation strongest, suggesting 
that only one of the two overlapping components strongly 
participates in the modulation cycle. 

At the pulse longitudes where the amplitude of the 
~ 20P modulation pattern is weakest the subpulse phase 
has the most complex pulse longitude dependence. In the 
first half of component II (pulse longitude ~ 86°) a rapid 
swing in subpulse phase occurs which cover almost a com- 
plete modulation cycle of ~ 20P (corresponding to 360° in 
subpulse phase). Similarly between components III and IV 
(around pulse longitude 107°) the modulation power disap- 
pears accompanied by a steep negative slope in the subpulse 
phase gradient. We note, however, that the power affected 
by this feature is very weak compared with the total power 
in the modulation pattern. 

The main conclusion of the subpulse phase analysis is 
that the regions where the ~ 20P modulation pattern is 
strongest a firm phase-locking is present between the MP 
and IP. We note that the delay accompanied by this phase- 
locking is found to be identical in completely different ob- 
servations recorded at different dates. This implies that the 
locking is sustained over periods much greater than the du- 
ration of a single observation and is likely to be a permanent 
feature of the pulsar. 



4 DISCUSSION 

4.1 The emission of PSR B1055-52 

4.1.1 Modulation features 

PSR B1055-52 is a pulsar with emission of great complexity. 
It is therefore important to assess the dominant features 
of the emission and their implications for possible physical 
models. 

Its most obvious feature is the ~ 20P modulation at 
both poles (Fig. [3} . Here we have demonstrated that this is 
not a precise periodicity, nor even a jitter about a specific 
periodicity, but loosely-connected quasi-periodic sequences. 
The sequences show no sign of 'drifting subpulses', and 
their intensity bands vary in intensity. In these respects 
PSR B1055-52 resembles PSR B1702-19 (jWeltevrede et all 
12007,) , also a near perpendicular rotator. So it may be argued 
that the lack of drift may be a consequence of being highly 
inclined. But regular mo dulations without drif t are also 
found in PSR B0656-hl6 (IWeltevrede et al.ll2006l), which is 
inclined at merely ~ 30° (|Everett fc WeisberdbOOll '). So the 
lack of drift may be simply a consequence of being relatively 
young (as many highly inclin ed pulsars are as according to 
IWehevrede fc Johnstonll2008l l . 

The second remarkable feature is the phase-locking of 
the IP and MP. Inspection of Fig. [5] shows that at times 
much (but not all) of the MP is imitating the IP in detail of 
intensity and pattern. The intrinsic delay of 3P is sufficiently 
large to be puzzling when measured in terms of light-travel 
distance (can the signal between the poles have travelled 
19 light-cylinder radii without losing coherence?). Yet the 
delay is small enough compared with the 20P modulation 
that we must conclude that information is conveyed from 
the IP to the MP, rather than the reverse. This would make 
the I P the 'driving' agent (as originally surmised by Biggsl 
Il990t ). and cross-correlations between its components show 
that despite its complex profile the IP effectively operates 
as a single component. This is illustrated by the subpulse 
phase plot of Fig. [T] where the IP modulation remains at a 
single phase, spread across all components. However, in the 
geometry established in Paper I, the IP is located on open 
polar field lines (as expected for pulsars' radio emission), 
and it is not clear how it can 'drive' the MP, which, under 
the assumption of a simple dipolar structure, is observed on 
putative closed field lines well within the light cylinder. 

The full picture represented by Fig. [7] is complex and 
potentially revealing. It is possible that the open field lines 
on the pole opposite to that of the IP may share the IP's 
modulation but these are anyway invisible to our sightline. 
Note that the MP, for reasons unknown, has a 'missing' lead- 
ing half as shown in Fig. 7 of Paper I. This implies that the 
closest we get to viewing the open field lines associated 
with the MP pole are the field lines on which the MP's 
leading components (I and II) lie, and with increasing pulse 
longitude our sightline moves progressively deeper into the 
'auroral' zone. Therefore component IV seems to define a 
boundary with a much-reduced corotating zone (see Fig. 8 of 
Paper I), although we should note that the implied presence 
of current flow within such a zone would require a modifi- 
cation of the assumed dipolar field lines. The ~ 20P modu- 
lation of the MP can be argued to appear at three different 
phases at three different longitudes. In components I and 
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the first part of component II the 20P modulation appears 
weakly and intermittently at zero lag with respect to that 
in the IP, close to the (hidden) open field lines of the second 
pole and therefore suggestive of synchronised two-pole mod- 
ulation. Moving deeper into the magnetosphere, and partly 
overlapping the earlier components, the second half of com- 
ponent II and component III are dominated by the princi- 
pal modulation at a lag of 2.5P (which corresponds to an 
intrinsic delay of 3P). A third phase of the 20P modulation 
appears in component IV, where the emission is actually al- 
most anti-correlated with the IP emission and appears in 
Fig. [7] to precede exact anti-correlation by ~ 40°. This effect 
can be seen, for example, around pulse 320 in the pulse- 
stacks of Fig. (2] Note that at two points in the profile (one 
just 2° before the pulse longitude of peak emission and one at 
pulse longitude ~ 106°) transitions from correlated to anti- 
correlated phases are marked by a rapid swing in subpulse 
phase. In any case, it is clear that the entire magnetosphere 
(with the possible exception of the small corotating region) 
is participating in the 20P modulation at multiple phases. 

It should be noted that this picture is not necessar- 
ily inconsistent with the two-pole model for PSR B1702-19 
j Weltevrede et al.ll2007l ). In that pulsar it was argued that 
modulated emission was intrinsically in phase at both poles 
and the emission was seen on open field lines from both 
poles. In PSR B1055-52 open field lines are only visible to 
us at the IP so that it is quite possible that the pole oppos- 
ing the IP is in exact phase with the IP. At the MP of PSR 
B1055-52 our sightline then catches a delayed feedback on 
field lines which are invisible to us in PSR B1702-19. 

There is also an interesting parallel with the behaviour 
of PSR B1822-09. Thi s pulsar is also mos t probably a near 
perpendicular rotator (|Backus et al.ll2010|') and has two dis- 
tinct modes of emission ('B' and 'Q': lFowler eTallligSlI ). In 
the B-mode the IP disappears and the MP gains an addi- 
tional component, so these component s are anti-correlated 
l|Fowler fc Wrightlll982l : [Gil et al.lll994l '). Even within the Q- 
mode the IP is anti-corr elated with one of t he MP subcom- 
ponents (see Fig. 11 of iBackus et al.ll2010t ). This is analo- 
gous to the behaviour of the IP and component IV in PSR 
B1055-52 and may have a common physical explanation. 

4.1.2 Pulse characteristics 

Figures [5] and E] give us very different views of the emis- 
sion. In the former we see a contrast between the spiky MP 
and broader IP emission, while in the latter we note a near- 
identical MP and IP energy distribution. The energy distri- 
bution's lognormal structur es (a common feature of pulsars; 
iBurke-Spolaor et al.l l|2012h ') imply a random multiplicative 
process, and the fact that they have the same standard de- 
viation appears to be telling us that identical processes are 
at work throughout the magnetosphere. 

However, we found no correlation between the MP and 
IP fiuxes at timescales longer than the modulation cycle, so 
the fact that the distributions have the same mean intensity 
(corresponding to the same MP and IP flux in the integrated 
profile) must be attributed to chance, and cannot favour our 
sightline and observing frequency and cannot be the result 
of reflections or refractions of emission from a single source. 
On the other hand, in the pulse-stacks of figure [2] we see the 
consequence of our particular sightline. So the appearance of 



a pulse, spiky or smeared, could depend on how our sightline 
samples the pulsars field lines, while the form of the energy 
distributions are related to the underlying processes. 

We also note that PSR B 1055-52 is not the only pulsar 
which shows radi o pulses with different characteristics in 
their MP and IP. iHankins fc Eilekl (|2007l ) have shown that 
the giant pulses observed in the IP of the Crab pulsar are 
much broader compared to those in the MP, reminiscent 
of the broader subpulses seen in the IP of PSR B1055-52. 
However, the IP pulses observed in the Crab pulsar show 
a very different observing frequency structure compared to 
those in the MP, while our data shows no such effect. Hence 
it is unclear if these two pulsars are linked or not. 



4.2 Seeking a physical model 

The essentially geometric results of Paper I, and those pre- 
sented hitherto in this paper, make no physical assumptions 
other than that we are observing a rotating dipole emitting 
polarized radiation parallel to its field lines. To develop a 
physical understanding, we need to explain (1) the origin 
of the dominant ~ 20P modulation cycle, (2) why there is 
a considerable phase-locked delay between the IP emission 
and that at the MP and (3) why this pulsar emits radiation 
on field lines which are normally deemed to be closed. 



4.2.1 Intrinsic models 

The polar cap in odel (|Ruderman fc Sutherland! Il975l : 
iGfl fc Sendvi3l2OO0l '). based on a dischar gmg vacuum gap 
above the poles, would have difficulty in accounting for the 
diverse emission of PSR B1055-52. There is no evidence ei- 
ther in the integrated profile or the single pulse behaviour to 
suggest an underlying conal structure and associated drift- 



ing sub pulses such as is found in so ma,ny, ge r ierally older, 
pulsars l|Rankinlll98i : IWeltevrede et al.ll2006l |2007| ). How- 
ever, it has been argued ( Gil fc Sendvkl 120001 ) that in rel- 



atively young pulsars the polar sparks will be narrow and 
spiky and display intensity modulation rather than drift. 
Spiky emission is indeed found at the MP of PSR B1055-52 
(though not at the IP). Nevertheless, all versions the po- 
lar cap model require the presence of non-dipolar field line 
structures close to the star's surface, possibly with a local 
sunspot-like structure, t o give sufficient field line curvature 
to trigger pair-creation (|Gil fc SendvkllioOOl ). However, this 
argument can be turned on its head: why would two very 
different local poles then have magnetic structures which 
generate precisely the same P3 periodicities? Furthermore, 
the poles share the same value of P3, even though its precise 
value varies intermittently around 20P. This, and the exact 
synchronised phase-lock delay of 2.5P, demands some form 
of inter-pole communication which the polar cap model does 
not describe in its present form. 

By contrast, a model of non-radial neu tron star oscilla- 
tions (adapted bv lClemens fc Rosen|[20o3 from the astero- 
seismology of magnetic white dwarfs) provides a natural way 
of inter-pole phase-locking via the surface of the neutron 
star. This has alre ady been pointed out in the similar case 
of PSR B1702-19 IWeltevrede et all 120071 '). In this model, 
the neutron surface is seen as having multiple narrow con- 
centric strips centred on the magnetic poles which oscillate 
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at a rapid frequency, with each at 180° phase shift from its 
neighbour. The 20P emission modulation then results from 
beats between the oscillation frequency and the observer's 
sampling rate of 1/P, and could be expected to be the same 
at both poles. The phase shift of the modulation is harder 
to explain: this should either be zero (as it is in the case 
of PSR B1702-19) or 180°, at least in the simplest version 
of the model. Neither a zero or 180° phase shift is consis- 
tent with the intrinsic 3P delay observed in PSR B1055-52. 
However, an advantage of this model is that it can plausi- 
bly allow for oscillations beyond the confines of the polar 
cap. The polar cap boundary is conventionally defined by 
the size of the light cylinder, and not by the neutron star's 
surface properties which would determine the oscillations. 
Thus emission in the closed magnetosphere is theoretically 
possible - but it would still be necessary to explain why PSR 
B105 5-52 is appa r ently so exceptional in this respect. Note 
that iKeith et al.] (|2010l ') have argued for emission outside 
the conventional open field line region for other interpulse 
pulsars. 

If we insist that information between the poles is 
magnetosphere-borne, the intrinsic delay time of ~ 3P is ex- 
traor dinarily long. In his original study of this pulsar, iBiggd 
I1990I pointed out that the light travel distance along a field 
line from one pole to the other is unlikely to be more than 
~ I.IP. An intrinsic delay time of ~ 3P corresponds to a 
light travel distance of about 19 light-cylinder radii. One way 
to achieve this would be for the polar emission to undergo 
refiections within the magnetosphere and/or some kind of 
magnetospher e- wide drift before becoming observab l e at the 
oppos ite pole ( Wright 2003; Me lrose fc Yu en '201 iV [Wrightl 
l|2003f ) invokes interactions between the pulsar's polar region 
and the null surfaces. In a highly-inclined pulsar the null sur- 
face can str etch to great heig hts above the pole along the ro- 
tation axis omani fc Wal ters 2010) and drifting particles 
may encounter electric gaps some 10 light-cylinder from the 
poles. Pair-generated particles would take a round trip time 
of the order of the intrinsic 3P delay time observed in PSR 
B 1055-52 and may well return to re gions w ell away from the 
polar gaps. In a more general paper iMelros c & Yuen (2011) 
argue strongly that corotating 'dead' zones cannot exist in 
an inclined rotating pulsar magnetosphere due to the in- 
nate inability of the pulsar to screen the inductive electric 
field. Hence particle drift must occur throughout the mag- 
netosphere and not only on open field lines. Although these 
authors only detail their model for nearly-aligned pulsars 
it has obvious wider application to PSR B1055-52 and its 
ideas would be well worth pursuing. 

4-2.2 Extrinsic models 

It seems that any model based on an isolated neutron star 
has difficulty in replicating the full complexity of the ob- 
servations. Could it be that the complexity arises from the 
intrusion of ex t raneo us material into the magnetosphere? 
ILuo fc Melros3 l|200'i1 ) have outlined such a model, i n which 
neutral material leaks into the corotation region (|Tsvganl 
Il977i : .Wright. 1979 ') and, following ionisation, modulates the 
radio emission through (ir)regular instabilities. Auroral cur- 
rents would cause the field lines to drift compared with coro- 
tation and thereby generate periodic modulations and quite 
plausibly generate a profile of multiple components. How- 



ever, one would expect both poles to be modulated at al- 
most the same phase (i.e much less than a single period), 
and, again, the rigid phase difference of 2.5P between poles 
remains hard to explain. 

A de bris disk akin to the i ntrusi on model has been sug- 
gested by ICordes fc ShannonI (|2008l ) as an explanation for 
intermittent m odulations of pulsar ra dio emission. Such a 
'fallback' disk (|Michel fc Desslerl|l98l| ) is argued to arise as 
neutral supernova material which has remained gravitation- 
ally bound to the neutron star gradually spirals inwards, 
spinning ever more rapidly as the star is approached, finally 
crossing the light cylinder. Individual rocks of various sizes 
are held together through their tensile strength and thus 
survive tidal disruption, despite gradual erosion through sur- 
face ionisation. For our present purposes this model has the 
advantage of introducing a natural source of periodic modu- 
lations at roughly the timescales required. In PSR B1055-52 
at just under half a light cylinder radius the orbital period 
of such material would be around 20 times the stars rotation 
period. Thus well within the light cylinder, yet outside the 
deduced closed region and perhaps sufficiently far from the 
neutron star to escape destruction by evaporation, we could 
expect the range of periods we observe as pulse modulation. 
The disk could comprise several distinct rings of fragmented 
neutral material (weak modulation), and/or be dominated 
by a small number of larger asteroids (strong modulation). 
An important advantage of this model is that an orbiting 
body might be viewed with an orbital time-delay at oppos- 
ing poles and thereby explain the phase-locked delay. Fur- 
thermore, as might be expected from a passing source, the 
~ 20P modulation appears in the pulse-stack (Fig. [Sjl as 
short-lived streaks of emission. 

The attractive features of the disk model made us fur- 
ther investigate its consequences. One critical expectation of 
the model is that the ~ 20P modulations should ultimately 
refiect an underlying rigid orbital periodicity, however much 
random jitter becomes imposed on it by the pulsar's mag- 
netosphere. We therefore took the brightest sequence of our 
data (the pulse-stack shown in Fig. and tested it system- 
atically with folds of the IP emission at different periodici- 
ties in the range 18 — 20P. As discussed in Section 13.31 the 
possibility of a fixed orbital timescale is ruled out. 



5 SUMMARY AND CONCLUSIONS 

In Paper I we were able to establish that PSR B1055-52, 
with its unusual and highly asymmetric integrated radio pro- 
file (Fig. [ij , belongs to the relatively small class of pulsars 
where we are viewing radio emission on field lines emanat- 
ing from both magnetic poles. Such pulsars present us with 
a rare opportunity to test emission theories, especially if, as 
here, their single-pulse emission exhibits specific modulation 
patterns at each pole. 

The radio emission of PSR B 1055-52 defies any tradi- 
tional pulsar category: there is little support for the picture 
of symmetrically-placed nested emission cones and an associ- 
ated pattern of subpulse drift. Instead we find a phase-locked 
delay between the modulations at opposite poles - and peri- 
odic modulation occurring on field lines normally attributed 
to a corotating 'dead' zone. Furthermore, the emission of vir- 
tually every component of both the MP and IP has a differ- 
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ent characteristic, ranging from smooth periodic to 'spiky' 
unmodulated. 

These results represent a major conundrum for pulsar 
theorists. The ~ 20P modulation found at both poles, al- 
though not precisely periodic, is phase-locked over the many 
epochs spanning our observations and therefore almost cer- 
tainly a permanent lock. In the profile components at both 
poles where this modulation is strongest, the pattern takes 
the form of one or two strong pulses followed by little emis- 
sion until the next band. In thi s respect, the emission re- 
sembles that of PSR B1702-19 () Weltevrede et al.ll2007l 'l. a 
pulsar similar in age and inclination - but whose interpulse 
peak separation is close to 180° and we can be confident we 
are observing emission on open field lines above the polar 
caps at alternate poles. By contrast, at one of the poles of 
PSR B1055-52 (the MP) we are observing emission on field 
lines which are well away from the polar cap and deep within 
what is normally deemed to be the closed co-rotating zone. 

The most striking feature of our present study is the 
magnitude of the phase-locked delay between the IP and 
MP modulations. At ~ 3P it implies relativistic commu- 
nication distances far exceeding the pulsar's light cylinder 
radius (corresponding to a light travel time of 0.16P), yet 
is short compared with the ~ 20P modulation it deliv- 
ers. In the case of PSR B1702-19 the intrinsic phase lock 
was virtually simultane ous and it was possible to argue 
l|Weltevrede et al.1 12007| ) that the poles may communicate 
to one another through the body of the neutron star itself 
(possibly via non-radial oscillations). Now this argument is 
much harder to make. 

Differences in the chance orientation of the observer's 
geometry may be responsible for the different appearances 
of PSRs B1702-19 and B1055-52. It is possible that both 
pulsars emit their quasi-periodic signals simultaneously at 
both poles on open field lines but are only observable in PSR 
B1702-19. These then later (in both pulsars) generate a re- 
sponse signal from certain closed field lines but are only ob- 
servable in PSR B1055-52. The simple uniform-phase nature 
of the IP's modulation compared to the complex multi-phase 
response of the MP (Fig. [7)l combined with the fact that the 
emission pattern of the MP is lagging that of the IP strongly 
suggests that the IP somehow drives the MP. We argue from 
our observations that the 20P bursts are magnetosphere- 
wide phenomena, synchronised at both poles and simulta- 
neously affecting regions well away from the poles. 

As an explanation for the phase-locking we initially 
found the idea of a circumstellar disk very appealing. An 
'asteroid' at a radius well within the light-cylinder can in- 
deed have an orbital period commensurate with the observed 
20P modulation, and this would be a natural way to bring 
about a phase-lock of significant delay. However, attempts 
to find the underlying orbital period by folding the modula- 
tion sequences met with failure: the ~ 20P modulations are 
not a simple jitter superimposed on precise periodicity. 

Th e popular (and in many cases convincing) polar cap 
model (|Ruderman fc Sutherlandlll975l : iGil fc Sendvkll20od ') 
has difficulties in explaining inter-pole locking. The model 
requires the emission modulation patterns to be fixed by lo- 
cal conditions at each magnetic pole and does not involve 
magnetosphere-wide interaction. Even if polar cap sparks 
can account for the emission behaviour of older 'drifting' 
pulsars, in the cases of PSRs B1055-52 and B1702-19 some- 



thing qui te different i s needed. Possibly a m odel akin to 
those of (| Wright! lioO^ : iMelrose fc Yuenll201ll 'l could work, 
linking separated regions of the magnetosphere and caus- 
ing the entire magnetosphere to participate in the subpulse 
drift. Alternatively, it is conceivable that irregularities in the 
pulsar's electrical current balance might be responsible. 

It should be noted that PSRs B1702-19 and B1055-52 
are not the only known examples of inter-pole communica- 
tion. The IP of PSR B18 22-09 was shown to int e ract with 
its M P many years ago (|Fowler fc Wrightl 1 19821 : iGil et al.1 
Il994 ) and its geometry has r ecently been c onfirm ed as a 
near-perpendicular rotator bv lBackus et al.l (|2010l ). In this 
pulsar a periodic modulation (~ 42P) is again shared by 
both the MP and IP, and the IP itself only appears dur- 
ing the one of th e pulsar's two modes. There is evidence 
(|Lvne et al.ll2010l ) that the mode-change is accompanied by 
a change in the spin-down rate. Clearly, the modulated radio 
emission of near-perpendicular rotating pulsars, although in- 
trinsically weak, is a diagnostic of powerful events. 

Taken together, the locked and subtle interactions be- 
tween the poles of PSRs B1055-52, B1702-19 and B1822-09 
must contain major clues to the nature of the elusive pulsar 
emission mechanism and the structure of the pulsar magne- 
tosphere. We have here focused on the single-pulse behaviour 
of PSR B1055-52 and briefiy considered a range of models. 
None satisfy all the observational constraints, so much work 
remains for those seeking to understand pulsars. 
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